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Abstract 

A theoretical  model  and  computations  of  the  I-V  characteristics  of  both  doped  and  undoped 
long  carbon  nanotubes  in  a strong  axial  dc-fields  at  room  temperatures  have  been  presented. 
Negative  differential  conductivity  has  been  predicted.  It  has  been  shown  that  \dI/dV\  for 
metal  carbon  nanotubes  in  the  region  of  the  negative  differential  conductivity  significantly 
exceeds  corresponding  values  for  semiconducting  ones.  The  predicted  effect  makes  possible 
the  design  of  wave-generating  nanotube-based  diods  for  submillimeter  and  infrared  ranges. 

1.  Introduction 

Since  the  discovery  by  Iijima  of  carbon  nanotubes  (CNs),  a great  deal  of  interest  has  been 
focused  on  these  quasi-one-dimensional  monomolecular  structures  because  of  their  unique  phys- 
ical properties  (mechanical,  electrical,  optical,  etc.)  and  the  rapid  experimental  progress  in  the 
controlled  preparation.  Processes  of  electron  transport  in  strong  external  fields  when  nonlinear 
effects  are  constitutive  are  of  a great  interest  for  potential  applications  in  nanoelectronics  and 
for  experimental  diagnostic  of  CN  themselves. 

The  current-voltage  (I-V)  characteristics  for  tunnelling  electrons  in  individual  single-wall 
CNs  at  law  temperatures  were  measured  in  Refs.[l,  2j.  At  such  temperatures  ksT  >C  £c,  A £ 
and  conduction  occurs  through  well  separated  discrete  electron  states;  here  fcg  is  the  Boltzmann 
constant,  T is  the  temperature,  £c  is  charging  energy,  A£  = ttHvf/L  is  the  energy  level  spacing 
with  vf  as  the  Fermi  speed  and  L as  the  CN  length.  At  the  above  condition  the  current 
is  produced  by  the  electrons  tunneling  through  CN  in  the  presence  of  the  Coulomb  blockade 
induced  by  the  long-ranged  (unscreened)  Coulomb  interaction. 

Current  instability  in  CNs  is  of  great  interest.  The  instability  appears  due  to  the  negative 
differential  conductivity  (NDC)  in  I-V  curves  of  CNs.  Nonlinear  coherent  transport  through 
doped  nanotube  junctions  was  considered  in  Ref.  [3].  It  was  also  shown  the  possibility  of  NDC  for 
tunelling  electrons.  In  Ref.[4]  a theoretical  phenomenological  analysis  of  the  I-V  characteristics 
of  undoped  CNs  at  room  temperatures,  when  kgT  > £c,  A£.  The  current  was  produced  by 
free  charge  carriers  (quasiparticles)  which  are  7r-electrons  moving  in  the  field  of  the  crystalline 
lattice.  The  nonlinearity  of  the  I-V  characteristic  appears  due  to  the  nonlinear  properties  of 
the  quasiparticle  gas  . The  negative  differential  conductivity  region  with  dl/dV  < 0 in  CNs  aws 
predicted  in  a certain  range  of  the  field  strength. 

2.  Theory 

This  report  extends  Ref.[4]  for  the  case  of  doped  CNs.  Let  us  consider  an  doped  single-wall 
zigzag  CN  (m,0)  exposed  to  a homogeneous  axial  dc-field  Ez,  Ez  = V/L,  where  V is  the  voltage 
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between  the  CN  ends.  We  shall  apply  the  semi-classical  approximation  considering  the  motion 
of  7r-electrons  as  a classical  motion  of  free  quasi-particles  with  dispersion  law  extracted  from  the 
quantum  theory  in  the  tight-binding  approximation  [5]. 

The  motion  of  quasi-particles  in  an  external  axial  electric  dc-field  is  described  by  the  Boltz- 
mann kinetic  equation 

e£,^  = -i[/(p)-F(p)l,  (1) 

where  e is  the  electron  charge,  F( p)  is  the  equilibrium  Fermi  distribution  function  and  r is 
the  relaxation  time.  The  relaxation  term  of  the  equation  (1)  describes  the  electron-phonon 
scattering  [6,  7],  electron-electron  collisions  etc. 

Utilizing  the  method  originally  developed  in  the  theory  of  quantum  semiconductor  super- 
lattices [8]  we  can  construct  an  exact  solution  of  kinetic  equation  (1)  without  assinning  the 
electric  field  to  be  weak.  First,  note  that  the  distribution  function  /( p)  is  periodic  in  pz  with 
period  2n/a,  a = 3b/2h,  b = 1,42A.  Then,  taking  into  account  the  transverse  quantization,  the 
distribution  function  can  be  presented  by: 

m 

/( P)  = &P<>>  S s(P<t>  ~ sAP4>)  £ frsemrPz , (2) 

s=l  r^O 


where  frs  are  coefficients  to  be  found,  S(x)  is  the  Dirac  delta- function,  A — 7t\/3 /am.  The 
equilibrium  distribution  function  F( p)  can  be  expanded  in  the  analogous  series  with  the  coeffi- 
cients as  follows 


Fr,  = 


a 
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2n/a . 
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1 + exp{(£s(pz)  - p)/kBT}  1 + exp{(-£s(p*)  - p)/kBT} 


dpz-  (3) 


In  this  equation  first  term  descirebes  the  contribution  of  conduction  band  and  the  second  one 
of  valence  band.  Here  p is  a chemical  potential  of  CN.  It  describes  doping  (p~0  for  undoped 
CNs  [4]).  p can  be  varied  within  wide  range  of  values.  Accordingly  to  [3]  p ~ 0.3  eV  for 
dopping  by  Au  substrate,  -0.5  eV  < p < 0.5  eV  for  doping  by  typical  alcali  and  halogen  atoms. 
For  the  KCs  doped  CNs  p ~ 2.0  eV  [9].  Substitution  of  both  expansions  into  Eq.  (1)  gives 
frs  = Frs/(  1 + irrO),  where  O = aeEz  is  the  Stark  frequency. 

The  surface  current  density  is  defined  by 


• ( p \ _ _ 8e7o  ^ r2ftr 
z)  y/Shbm  1 + (rfir )2 


m 

F rs^rs- 

3 = 1 


(4) 


Here  70  ~ 3.0eV  is  the  overlapping  integral.  This  equation  states  the  basis  for  the  evaluation  of  I- 
V characteristics.  As  it  has  been  stated  above,  the  quasi-particles  motion  is  described  classically. 
Thus,  both  interband  transitions  and  quantum-mechanical  corrections  to  the  intraband  motion 
are  left  out  of  account  in  this  model.  This  imposes  the  limitation  on  the  external  electric  field 
strength:  \EZ\  < 70/2 eR. 

The  Coulomb  electron-electron  interaction  has  been  also  left  out  of  account  in  our  approach. 
The  role  of  this  mechanism  as  applied  to  CNs  was  considered  in  a number  of  papers,  see  [10]  for 
example.  It  has  been  established  that  the  short-range  electron-electron  interaction,  typical  for 
CN  arrays  (’ropes’),  have  only  a weak  effects  at  high  temperatures. 


3.  Numerical  Results 

The  I-V  characteristics  obtained  via  numerical  calculation  of  Eq.(4)  are  presented  in  Fig.l  and 
Fig.  2 for  metal  (m  = 3 <7,  q is  an  integer)  and  Fig.  3 and  Fig.  4 for  semiconducting  (m  7^  3 q) 
zigzag  CNs. 


jz1  o',  A/cm 
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Fig.  1 Undoped  metall  zigzag  CNs. 
T=287,5°K,  t=3*1  O'12. 


Fig.  2 Doped  zigzag  CNs  with  m=3q  (q  is  an  integer). 
H=0.2  eV,  T=287,5°K,  -c=3*10‘12 


The  figures  show  the  linear  dependence  of  jz  on  Ez  at  weak  strengths  of  the  external  field 
both  for  doped  and  undoped  CNs;  it  corresponds  to  the  region  of  ohmic  conductivity.  As  Ez 
increases,  the  value  djz/dEz  growth  smaller  and  at  Ez  = the  current  density  reaches  the 

maximum  value  j™ax.  Further  increase  of  Ez  results  in  the  decrease  of  jz.  Thus,  predicted  in  [4] 
the  region  with  the  negative  differential  conductivity  ( djz/dEz  < 0),  in  the  I-V  characteristics 
of  undoped  CNs  can  be  also  observed  in  the  case  of  doped  CNs. 

Fig.l  and  Fig.  3 demonstrate  that  E^nax'>  depends  on  neither  number  m nor  the  conductivity 
type  (metal  or  semiconductor),  whereas,  j™ax  shows  the  different  dependencies  on  m for  undoped 
metal  and  semiconducting  CNs.  For  metal  CNs,  j™ax  decreases  with  m while  it  increases  for 
semiconducting  ones.  As  m ->■  oo,  j™ax  for  metal  and  semiconducting  CNs  tends  to  the  same 
limit  from  opposite  sides.  Generally,  at  large  m,  the  I-V  characteristics  of  different  CNs  are 
coming  close  and  in  the  limit  case  m ->•  oo  they  reduce  to  I-V  characteristic  of  the  plane  graphite 
monolayer.  It  should  be  noted  that  the  metal  CNs  exhibit  much  lager  NDC  as  compared  to 
semiconducting  ones.  Fig.  2 shows  that  the  doping  of  metal  CNs  makes  J%iax  some  times 
larger.  Doping  of  semiconducting  CNs  with  small  m does  not  lead  to  noticeable  changes  but  as 
m increses,  current  through  doped  doped  semiconducting  CNs  increases  and  values  of  current 
density  become  one  order  with  ones  for  conducting  undoped  CNs  (Fig.  4). 

4.  Conclusion 

In  Summary,  we  have  predicted  the  NDC  effect  in  both  doped  and  undoped  CNs,  which  is 
expected  to  be  observable  in  sufficiently  long  CNs  at  room  temperatures.  I-V  curves  are  expected 
to  be  effectevely  controlled  due  to  doping. 
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Fig.  3.  Undopped  semiconducting  zigzag  CNs. 
T=287,5°K,  t=3*1 0"12. 
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Fig.  4.  Doped  zigzag  CNs  with  m=3q+1  (q  is  an  integer). 
H=0.2  eV,  T=287,5°K,  -c=3*10'12 
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